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In order to investigate charge ordering phenomena under electric field, static nonequilib- 
rium Hartree approximation (SNHA) method is formulated on the basis of the nonequilibrium 
Green's functions introduced by Keldysh. By applying the SNHA to the 3/4- filling extended 
Hubbard model on anisotropic triangular lattice, we study the stabilities and amplitudes of 
3-fold and horizontal charge orders in and ^cZ-(BEDT-TTF)2X salts under the electric field. 
The obtained results show that the electric field stabilizes the 3-fold state in comparison to the 
horizontal state. The amplitude of the 3-fold state tends to decrease by the field, whereas that 
of the horizontal state does not change. 

KEYWORDS: charge order, anisotropic triangular lattice, organic conductor, static nonequilibrium condi- 
tion, Hartree approximation, nonequilibrium Green's function 



One of the families of organic crystals, (BEDT- 
TTF)2X (in short (ET)2X), has quasi-two- dimensional 
nature. Its donors, BEDT-TTF (ET) molecules, and 
acceptors, Xs, form monolayers respectively, and they 
are stacked layer by layer. Spatial arrangements of ET 
molecules can be controlled by the anion molecules and 
hydrostatic pressure.^ Among the different polytypes, 
ET layers in and Od-^JV^ crystals, which are focused 
on in this paper, form 3/4-filling anisotropic triangular 
lattices: the former is metallic and has no dimerization, 
whereas the latter is insulating and dimerized. 

When the acceptor, X, is RbM'(SCN)4 {M' =Co, Zn), 
(ET)2X is a metallic ^-type material at room tempera- 
ture. When this Rb salt is cooled slowly (slow-cooling 
condition), it undergoes a structural phase transition 
from 9 type to Od at about 190K(= Tmi), and be- 
comes insulator. Accompanying this transition, the 
electronic state also changes: above Tmi, a metallic short- 
range charge order (CO) characterized by a diffuse rod, 
q = (1/3, /c, 1/4), is observed in X-ray measurements,^ 
while below Tmi, an insulating long-range CO character- 
ized by a diffuse rod, q = (0, /c, 1/2) = qi, is observed in 
NMR, X-ray, Raman scattering, and optical conductivity 
experiments. The latter CO is called the horizontal 
CO, which has two- fold periodicity in real space. 

On the other hand, when X is CsM'(SCN)4 {M' =Co, 
Zn), the crystal remains 9 type down to low tempera- 
ture and no phase transition is observed.^ Above 20K, 
it is metallic and a short-range 3-fold diffuse rod, q = 
(2/3, A:, 1/3) = q2, is observed. At about T = 20K, the 
resistance starts to increase by 10"^"^ times accompa- 
nied by an appearance of the short-range diffuse rod at 
q = qi, which coexists with q2.^^'^^ The wave number, 
qi, is same as that observed in Rb salts below Tmi i.e., 
the horizontal CO. Since the horizontal CO will be ac- 
companied by the Od type lattice structure, Cs salts are 
supposed to be in a super-cooled glass state consisting 
of inhomogeneous mixture of 9 and Od-^yv^ lattice struc- 
tures below 20K.^'i^ 
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Recently, surprising current- volt age {I-V) characteris- 
tics are observed in Cs salts below 20K, where the short- 
range 3-fold and the horizontal CO coexist: a gigantic 
nonlinear conductivity^^' and a thyristor like I-V char- 
acteristic.^^ Furthermore, the X-ray measurements under 
a constant current below 5K show that the intensity of 
the horizontal diffuse rod decreases whereas that of the 
3-fold remains. Therefore, the observed anomalous I-V 
characteristics are assumed to be due to the different de- 
pendences of CO state on the electric field or current. 
In Rb salts, under a rapid-cooling condition, the same 
phenomena as Cs salts have been observed. 

The main purpose of this paper is to understand these 
phenomena using an extended Hubbard model which 
is a canonical model describing charge ordered organic 
conductors. The preceding theoretical researches about 
charge ordering phenomena are often focused on stability 
and phase diagram under equilibrium condition. ^^^^ In 
this paper, we formulate static nonequilibrium Hartree 
approximation (SNHA) using Green's functions intro- 
duced by Keldysh. In this method, the free energy in 
the presence of static electric field can be calculated. The 
I-V characteristic is also investigated. We will show that 
the electric field stabilizes the 3-fold CO when the hori- 
zontal and 3-fold states are quasi degenerate. This may 
explain the X-ray measurements under a constant cur- 
rent. 

The model for the calculations is based on the fol- 
lowing three hypotheses. (1) All calculations are carried 
out in the or ^^^-type crystal, although Cs salts may 
have a glass state consisting of the 3-fold state with the 
^-type lattice and the horizontal state with the ^^^-type 
lattice structure. (2) The metallic state under the elec- 
tric field is assumed to be the 3-fold state in the ^c^-type 
crystal in order to compare the energies, although ex- 
perimentally, the crystal structure transforms into a new 
structure under the electric field, which is slightly differ- 
ent both from the and 6>^i-type. (3) We compare the 
energies of the 3-fold and horizontal state calculated in 
the uniform states. This means that the boundary ener- 
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Fig. 1. Crystal structure and transfer integrals of (a) 6 and (b) 
^(i-(ET)2X in a donor layer. 



Table I. Transfer integrals of 6>d-(ET)2RbZn(SCN)4 at 90K cal- 
culated by the extended Hiickel method.^ 



Transfer integral 
(xlO-^eV) 



cl 



6'd-(ET)2RbZn(SCN)4 
c2 pi p2 p3 



p4 



1.5 5.2 16.9 -6.5 2.2 



-12.3 



gies between the metallic and the insulating regions are 
neglected. 

Our Hamiltonian is written as the following extended 
Hubbard model: 



(1) 



where txy is the transfer integral between sites, x and 
and U and Vxy are on-site and nearest-neighbor Coulomb 
repulsion, respectively. The summation, ..\ , is 
over the pairs of nearest-neighbor sites. The operators, 
{cxcj) and Uxa represent a creation (an annihilation) 
and a number operator at site, x, with spin, a. Through- 
out our calculations, U and the ratio Vc/U^ are fixed at 
U = 0.7eV and Vc/U = 0.3, respectively, which are typi- 
cal values for organic conductors. The transfer integrals, 
txy^ are tc = 4.3meV and tp = ll.SmeV for the ^-type 
crystal as shown in Fig. 1. (a), and those for the ^c^-type 
crystal are given in Fig. 1. (b) and Table I. 

When we apply the Hartree approximation, the Hamil- 
tonian, is transformed into two parts, Hmf and a 
constant term. The former is bilinear in the electron op- 
erators. Taking twelve molecular sites as a unit cell, Hmf 
can be diagonalized in momentum space as, 



(2) 



kacT 



Here a new annihilation (creation) operator is written as 
Ckacr = X]/3 ^k/3Q;^k/5cr7 whcrc a and (3 represent the band 
index and the site number in a unit cell, and iXk/3a and 
£ka are the eigenvector and the eigenvalue of i^MF, re- 
spectively. In this paper, a unit cell includes 12 molecules 
in order that both 3-fold and horizontal CO are fit in the 
unit cell. On the other hand, when we study the ^-type 
crystal, a unit cell with 3 molecules is used, because the 
^-type crystal has no solution of the horizontal CO in 
the realistic parameter region. 

In the SNHA method that we formulate, observables 
and self-consistent equations can be derived in terms of 



the Green's functions introduced by Keldysh.^^ Here, we 
use the following Green's functions, so-called the lesser 
and greater Green's functions in the real space: 



G<y,{T,T)=i{cl,{t')c.a{t)) 



(3) 



with T = t — t' and T = ^(t + 1'), where the dependence 
of G on T can be ignored. In these definitions, Cxcr{t) 
(cj.^(t)) is the annihilation (creation) operator of Hmf, 
which also includes the electric field term, — ena^(^)x • E, 
and the impurity scattering term, both of which are bi- 
linear. In terms of the lesser Green's function in Eq. (3), 
(Hmf) can be obtained as. 



(Hmf) 



lim lim 



E 

xa 



dr "'y" 



-I lim > 

kaa 



r]- 



/duo 
2^' 



(4) 



o;e-''G<,(k,a;), 



where G^^(k, cj) is the Fourier transform of G^y^{r). 

Then, according to the theory of quantum kinetic 
equation based on the Keldysh's method, the lesser 
Green's function, G^^^^(k, lo) in the first order of |E| can 
be obtained as: 

G<^\k,u;) -m(k,a;){nF(a;) - Ak<(a;)}, (5) 

where np^uj) represents the Fermi distribution function. 
In (5), a(k, u;) and Ak are, 

7 



a(k, cj) = 



Ak = eE . 



m 7 



(6) 



(7) 



where the effective electron mass, m, in (ET)2X is as- 
sumed to be the free electron mass. In this formulation,^'' 
the relaxation rate, 7, naturally appears in a(k, cj) and 
Ak due to the presence of impurity scattering. In our cal- 
culation, the relaxation time, r = h/'j is fixed at 10~^s 
in order to reproduce the resistivity. 

We have extended the calculation for G^*^^^ (k, lo) up to 
the second order with respect to |E|. After some straight- 
forward manipulations, we obtain, 

(8) 

which we approximate into the following expression: 



G^^(k, uj) ia(k, uj)nF{uJ - Ak). 



(9) 



Details of calculations will be published elsewhere. Since 
7 {r^ 10~ ''eV), is negligible compared to kBT{< 10~^eV), 
we approximate a(k, cj) with delta function, 27rS{(jj — 
^kacr)- With use of the lesser Green's function in (9), 
{Hmf) becomes. 



MF) 



kaa 



jn^iuj - Ak). 



(10) 



In the same way, the self-consistent equations and expec- 
tation value of free energy and current are obtained in 
terms of G^^(k, u;). 

Figure 2 shows the obtained ground state energies of 
horizontal and 3-fold CO in the 6d type crystal for the 
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Fig. 2. Electric field dependence of the ground state energies at 
Vp/Vc = 0.9. The 3-fold CO state is stabilized by the electric field 
and it has a lower energy than the horizontal state for |Ec| > 
3.2V/m 
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Fig. 4. Intensity, /, of a diffuse rod, q, refiects both intensity 
(upper panels) and area of charge modulation (lower panels). 
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Fig. 3. The phase diagram against the temperature and the elec- 
tric field. The small solid and open circles in each region represent 
the charge distribution of each CO. 



case with Vp/Vc = 0.9. If the electric field is applied more 
than the critical field, Ec, the 3-fold state is stabilized in 
comparison to the horizontal state. We also calculated 
the gournd state energies for 0.8 < Vp/Vc < 1-0. When 
Vp/Vc is larger, the 3-fold CO state is more stable than 
the horizontal CO state. 

Next, we calculate the free energies of two CO states 
at finite temperature, T, ranged from 5K to lOK. The 
anisotropy, Vp/Vc^ is fixed at 0.89998 so that the en- 
ergetic competition between the 3-fold and horizontal 
states is clearly seen. Figure 3 shows the critical elec- 
tric field for various temperatures, which gives a phase 
diagram. As T increases, the 3-fold state is stabilized. 
Note that the free energy comparison is carried out in the 
6d-type structure, although experimentally the metal- 
lic 3-fold state is in the ^-type structure. In order to 
check the differences between the 3-fold state in and 
Od-type structure, we have studied the free energy of the 
3-fold state in the 6>-type, and found that the electric- 
field dependences are similar quantitatively. Therefore, 
if we compare the energies between the 3-fold state in 
the 6>-type, we expect to obtain the same phase diagram 
as in Fig. 3. 

We have also investigated the amplitude of the hori- 
zontal and 3-fold CO as a function of the electric field 
(not shown here). In this case we studied that the CO in 
the ^-type crystal instead of Od type, because the 3-fold 
CO in the Od type has the horizontal modulation. The 



obtained results show that the electric field dependence 
of the amplitude of the 3-fold CO is / = — a^^ + l.O, with 
a = 1.03 X 10"^ at T = 9.5K. The more the electric field 
is applied, the weaker the intensity becomes. If the elec- 
tric field is about 500 V/m, where nonlinear conductivity 
is experimentally observed, we can expect 26% reduction 
of the amplitude, although our formulation can not be 
applied in such a strong field. On the other hand, the 
amplitude of the horizontal CO does not change by the 
electric field. 

These results suggest that the electric field stabilizes 
the 3-fold CO but decrease its amplitude. We expect that 
the experimentally observed intensity of the diffuse rod, 
q2, is affected by these two effects. Let us assume two 
cases as shown in Fig. 4. In (a), the amplitude of the 
3-fold CO is twice stronger than that in (b). However, 
its area is half of that of (b). If one measures the X-ray 
diffraction in these two cases, one may observe diffuse 
rods with roughly the same intensity. Although it is nec- 
essary to consider the domain size and the boundary en- 
ergies, it is possible that the total intensity of the 3-fold 
CO state does not change under the electric field, while 
that of the horizontal state decreases. If this is the case, 
the present calculation explains the tendency of what are 
observed in the X-ray measurements by Sawano et al.^^ 

Finally, we would discuss the validity of the SNHA 
method. In the region of |E| < 1.0 V/m, the shift Ak in 
the nonequilibrium Green's function in Eq. (7) is about 
O.Oltc, which is about 1/100 of the band gap of Od- 
(ET)2X obtained in the tight-binding approximation. Al- 
though Ak is small, there is an error in the replacement 
of G<(2)(k,cj) in (8) into G<{k,u;) in (9). This error can 
be estimated by. 



G<(2)(k,M)-G<(k,M) 



(11) 



G<(2)(k,M) 

which is about 10~^ for |E| < l.OV/m. However, when 
|E| is larger than l.OV/m, the error exceeds 0.3, since 
/3Ak appearing in nri^ — Ak) becomes close to 1 at 
the low temperature (5K, for example). This means that 
in the higher electric fields, the error becomes larger, 
and higher order correlations in (8) become necessary. 
This is the reason why we only applied the SNHA for 
|E| < 6.0V/m in Fig. 2 and Fig. 3. To deal with the 
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Electric Field [V/m] 

Fig. 5. I-V characteristics of 3-fold CO in and 0^ type crystals 
show linear conductivity. 



Since the SNHA method cannot be apphed to higher elec- 
tric field, we did not discuss the gigantic nonlinear con- 
ductivity and thylistor-like conductivity observed exper- 
imentally. Although further improvements in the SNHA 
method are necessary to deal with the high-field prob- 
lems, the present method predicts tendency of behaviors 
of CO states under the static electric field. 
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Fig. 6. Three parameter regions give linear conductivity, non- 
linear conductivity, and tylistor-like characteristic, respectively. 



higher electric fields, further improvements are necessary. 

When we calculate the I-V characteristics in the 
present formalism, we obtain the linear conductivity as 
shown in Fig. 5. This is reasonable because the region 
that we study corresponds to the region where the linear 
conductivity is observed experimentally. The parameter 
regions are shown in Fig. 6 schematically. The present 
calculation corresponds to the left-hand side of this di- 
agram. In this linear-conductivity region, the calculated 
resistivities of O{10~^)rtcm in the 6d-type crystal are 
consistent with the measured resistivities of Rb salts of 
BEDT-TTF at about 200K. The calculated resistivities 
in the 3-fold CO of O{10-^)ncm in the l9-type are 1/10 
of the measured resistivities of Cs salts at about 200K. 

In summary, we formulate the static nonequilib- 
rium Hartree approximation (SNHA) method based on 
Keldysh's Green's functions in order to calculate the 
observables under the constant electric field. By ap- 
plying this method to the 3/4-filled extended Hubbard 
model, we study the ground state energy with respect 
to anisotropy, Vp/Vc^ the free energy at various tem- 
peratures, charge disproportion in a unit cell, and I-V 
characteristics. The results show that the electric field 
stabilizes the 3-fold state in comparison to the horizon- 
tal state. The charge amplitude of the 3-fold CO is de- 
creased by the applied field, while the horizontal CO 
state is not affected. We speculate that these results ex- 
plain the X-ray measurement in 6>-(ET)2X at low temper- 
ature, in which the diffuse rod, qi, has tendency to dis- 
appear, whereas q2 remains. The I-V characteristics for 
|E| < 6.0 V/m is almost consistent with the experiments. 
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